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Abstract. The solar corona has been revealed in the past decade to be a highly dynamic nonequi- 
librium plasma environment. Both the loop-filled coronal base and the extended acceleration region 
of the solar wind appear to be strongly turbulent, but direct observational evidence for a cascade of 
fluctuation energy from large to small scales is lacking. In this paper I will review the observations 
of wavelike motions in the corona over a wide range of scales, as well as the macroscopic effects of 
wave-particle interactions such as preferential ion heating. I will also present a summary of recent 
theoretical modeling efforts that seem to explain the time-steady properties of the corona (and the 
fast and slow solar wind) in terms of an anisotropic MHD cascade driven by the partial reflection 
of low-frequency Alfven waves propagating along the superradially expanding solar magnetic field. 
Complete theoretical models are difficult to construct, though, because many of the proposed phys- 
ical processes act on a multiplicity of spatial scales (from centimeters to solar radii) with feedback 
effects not yet well understood. This paper is thus a progress report on various attempts to couple 
these disparate scales. 
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INTRODUCTION 

Astronomers often define "turbulence" rather loosely, as either a collection of motions 
that are unresolved either spatially or temporally, or as dynamical oscillations that 
exhibit a broad-band spectrum of frequencies and no clear dominant frequency. Our 
knowledge about turbulence in the solar corona comes mainly from observations of this 
kind. Iron-clad evidence for the existence of a turbulent cascade in the corona awaits 
in situ exploration such as a "Solar Probe" could provide [|32|l . Substantial progress has 
been made, though, on the basis of remote- sensing observations, extrapolation inward 
from existing in situ measurements, and theoretical modeling. This paper summarizes 
a cross section of these recent efforts to better understand the role of turbulence in the 
corona. 

OBSERVATIONAL EVIDENCE 



The turbulent solar photosphere is the natural lower boundary condition for fluctua- 
tions higher in the atmosphere [37]. The photosphere displays a superposition of both 
quasi-laminar granulation (i.e., overturning convective cells) and smaller-scale stochas- 
tic motions in the dark intergranular lanes. The latter are also associated with 100 km 
sized concentrations of magnetic field (~1.5 kG) — known as G-band bright points or 
magnetic bright points — that are shaken transversely by the granulation and appear to 
contain enough energy to give rise to coronal Alfven waves [7J. 



Higher in the corona, plasma fluctuations reveal themselves by presenting variations 
in density, velocity, and magnetic field strength. These remotely measured quantities 
are most sensitive to the wave modes that carry the most energy, which are generally 
dominated by the longest wavelengths. Thus, direct measurements are typically inter- 
preted as ideal MHD fluctuations. (Some indirect measurements provide constraints on 
small-scale kinetic modes; see below.) The major detection techniques are listed here: 

1. Intensity modulations mainly probe variations in density, with the observed fluc- 
tuations being proportional to either p or (integrated along the line-of-sight) 
depending on the spectral band or lines used. Intensity oscillations measured with 
various instruments aboard the Solar and Heliospheric Observatory (SOHO) have 
implied the presence of compressive MHD waves channeled along magnetic flux 
tubes in coronal holes [9, 34]. 

2. Motion tracking in sequences of images can provide information about velocity 
fluctuations in the "plane of the sky" (i.e., the plane perpendicular to the line-of- 
sight direction). Cross -correlation techniques have been used to obtain the bulk 
solar wind acceleration of low-contrast "blobs" [42]. Recently, wavelet-enhanced 
prl] images and movies from the EIT and LASCO instruments on SOHO have 
enabled the filtering of fine-scale variations that were previously "hidden" in diffuse 
larger-scale coronal features. 

3. Doppler shifts generally allow velocities along the line-of-sight direction to be 
probed. Time-resolved sinusoidal oscillations have been seen in some coronal 



structures [|46L [350. Most often, though, waves that reach into the solar wind are 
diagnosed via "Doppler broadening" of spectral lines that arises because of averag- 
ing over the oscillating redshifts and blueshifts. In the low corona, where all plasma 
species are expected to be coUisionally coupled (and have identical temperatures), 
it is relatively straightforward to extract the "nonthermal" wave broadening from 
the thermal motions that also contribute to emission line widths. In the coUisionless 
extended corona, though, there appears to be at least a mild decoupling between Te 
and Tp as well as stronger preferential heating for heavy ions (see below). This 
complicates the analysis, but realistic limits can still be placed on transverse wave 



amplitudes [for summaries of existing data, see 11,0, [13, 



m. 



4. Radio sounding probes the conditions in the corona by measuring distortions in 
the net refractive index of plasma that intervenes between a receiver and either 
a spacecraft beacon or a cosmic source (e.g., a pulsar or radio galaxy). These 
distortions are sensitive to density (from scintillations), velocity (from drifting 
diffraction patterns), and the magnetic field (from Faraday rotation). The wide 
range of spatial scales probed by radio diagnostics have allowed new constraints 
to be placed on the properties of turbulence in the corona [2, 4, 2^1221], but there 



is still something of a disconnect between the quantities that are measured directly 
and the quantities predicted by theory. 

If the coronal fluctuations merely propagated upward on open field lines — without 
interacting with the background plasma — they would not be of much interest. Observa- 
tions of how the plasma is heated and accelerated, presumably by wave-particle inter- 
actions in the coUisionless outer corona, are useful as an indirect means of determining 



which wave modes are generated and damped. The Ultraviolet Coronagraph Spectrom- 
eter (UVCS) on SOHO measured extremely high temperatures of heavy ions, faster ion 
outflow compared to protons, and strong anisotropics (in the sense > 7||) for ion ve- 
locity distributions in the extended corona [ 25ll26ll27l]. These properties are shared by 
high-speed solar wind streams measured in situ [e.g.. 13011. 

The UVCS observations have led to a resurgence of interest in ion cyclotron reso- 
nance as a likely mechanism for producing this kind of preferential ion energization, 
and possibly also for heating the bulk plasma as well 10, 112, l27|] ■ Alfven waves that are 
ion cyclotron resonant in the corona have frequencies in the 10^-10^ Hz range, whereas 
the measured and inferred frequencies are much smaller; 10^^-10^^ Hz. Thus, turbulent 
cascade has long been considered a natural way to produce power at high frequencies 
from waves initially at lower frequencies. It is well known, though, that both numeri- 
cal simulations and analytic descriptions of MHD turbulence (with a strong background 
"guide field" like in the corona) indicate that the cascade from small to large wavenum- 
ber occurs most efficiently for modes that do not increase in frequency (i.e., primarily 
a fast cascade in and negligible transport in k\\ ~ (u/Va)- In the low-j8 corona, this 
type of quasi-two-dimensional cascade would lead to kinetic Alfven waves and prefer- 
ential electron heating (in 7j|), which is not observed. This issue remains a topic of active 
research, with several possible outcomes depending on the behavior of the anisotropic 
cascade when kinetic processes become important I13l[6l[22i. i29ll24l [39ll. 



CORONAL HEATING AND SOLAR WIND ACCELERATION 

Much of the above work dealt with determining the properties of the "microscopic" 
fluctuations that dissipate to heat the particles in the corona. Complementary progress 
has been made in constraining the "macroscopic" properties of the MHD turbulence that 
should eventually cascade down to the microscopic kinetic scales. In fact, if the precise 
distribution of energy into various channels (i.e., ^ Tp and 7]| 7^ T^) is ignored, it can 
be argued that the macroscopic level is all that is needed to compute how much energy 
will eventually be dissipated. This approach has been taken in recent models of turbulent 
heating of both closed loops in the low corona [38] and open flux tubes that reach into 
the solar wind [vLlsLlil. 

The remainder of this paper presents results from a self-consistent treatment of coro- 
nal heating and solar wind acceleration that used a phenomenological description of 
imbalanced MHD turbulence [jSt]. The only input parameters to these models were the 
photospheric lower boundary conditions for the wave spectra and the radial dependence 
of the background magnetic field along the flux tube. The models are the first self- 
consistent solutions that combine: (a) chromospheric heating driven by an empirically 
guided acoustic wave/shock spectrum, (b) coronal heating from Alfven waves that have 
been partially reflected, then damped via a turbulent cascade, and (c) solar wind acceler- 
ation from gradients of gas pressure, acoustic wave pressure, and Alfven wave pressure. 

The majority of heating in these models comes from the turbulent dissipation of 
partially reflected Alfven waves [see also|3ll[lli|45(]. Measurements of G-band bright 
points in the photosphere were used to specify the Alfven wave power spectrum at the 
lower boundary. Non-WKB wave transport equations were then solved to determine the 



degree of linear reflection due to radial gradients in the background plasma parameters 
(mainly the Alfven speed Va). The resulting values of the Elsasser amplitudes Z±, which 
denote the energy contained by upward (Z_) and downward (2+) propagating waves, 
were then used to constrain the energy flux in the cascade. We used a phenomenological 
form for the nonlinear transport that has evolved from studies of reduced MHD and 
comparisons with numerical simulations. The adopted volumetric heating rate (erg s^^ 
cm^^) is given by 

[e.g., 23, 50, 11, sj]. The transverse length scale represents an effective perpendicular 
correlation length of the turbulence, and we used the standard assumption that scales 
with the cross-sectional width of the flux tube [21]. The term in parentheses above is 
an efficiency factor that accounts for situations when the cascade does not have time 
to develop before the waves or the wind carry away the energy [10]. The classical 
Kolmogorov-like cascade is thus "quenched" when the nonlinear eddy time scale ?eddy 
becomes much longer than the macroscopic wave reflection time scale Ufi . In most of 
the models we used n = \ based on analytic and numerical models |T2, la], but we also 
tried n = 2 to explore a stronger form of this quenching. 

Figure 1 summarizes the results of varying the magnetic field properties while keeping 
the lower boundary conditions fixed [8] . The models included polar coronal holes, equa- 
torial streamers (as well as the full latitudinal variation between the two at solar mini- 
mum) and open flux tubes rooted in active regions. We found that a realistic variation 
of asymptotic solar wind conditions can be produced by varying only the background 
magnetic field geometry, as predicted by Wang and Sheeley [43, |48i |49] . Specifically, 
the models show general agreement with some well known empirical correlations: i.e., 
a larger coronal expansion factor gives rise to a slower and denser wind, less intense 
Alfvenic fluctuations at 1 AU, and larger values of both the 0^+/0^+ charge state ratio 
and the FlP-sensitive Fe/0 abundance ratio. Satisfying these kinds of observational scal- 
ings are necessary but not sufficient conditions for validating the idea that the solar wind 
is driven by a combination of MHD turbulence and non-WKB Alfven wave reflection. 

The models shown in Figure 1 are limited by being one-dimensional, time- 
independent, and one-fluid descriptions of a plasma that, in reality, is none of those 
things. Future work must involve including the divergent temperatures and flow speeds 
of protons, electrons, and various heavy ion species in the extended corona and helio- 
sphere. Even simple two-fluid effects (Je 7^ Tp) can affect the macroscopic distribution 
of heat flux [^, the dynamical stability of closed-field regions like streamers fil3i] . and 
possibly even the phenomenological form of the MHD turbulent cascade [[l6|]. Also, the 
inclusion of compressive MHD modes may also affect many important properties of the 
anisotropic cascade [e.g., 3]. 

It is also likely that in certain regions of the corona, waves and turbulence (by 
themselves) cannot be the whole story. For example, it seems increasingly clear that 
the bright coronal loops seen in UV and X-ray images are heated by some variety of 
intermittent magnetic reconnection. The question remains whether such a collection 
of bursty heating events is powered mainly by: (a) direct stressing of the magnetic 
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FIGURE 1. Summary of recent turbulence-driven coronal heating and solar wind acceleration models, 
(a) Adopted magnetic field geometry [T], with radii of wave-modified critical points marked by open 
diamonds, (b) Latitudinal dependence of outflow speed at ~2 AU for models with n = 1 (thick solid curve) 
and n = 2 (dashed curve), compared with data from the Ulysses polar pass in 1994-1995 (thin solid curve) 
ITtIi . (c) Radial temperature dependence of polar coronal hole model (solid curve), equatorial streamer- 
edge model (dashed curve), and strong-field active region model (dotted curve). Additional details about 
these models can be found in 



footpoints [e.g., uM, (b) newly emerging magnetic flux from below the photosphere 
[ 15ll40|] . (c) long-time buildup of non-potential shear [|28|] . or (d) some combination of 
the above ideas. It is possible that concepts from turbulence theory may be applicable to 
these kinds of processes as well [44, 33]. 
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